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Abstract 

We  present  the  results  of  semiempirical  quantum  chemical  calculations  on 
oligomers  of  poiy(p-pyridyl  vinylene)  (PPyV)  and  poly(p- pyridine)  (PPy). 

The  presence  of  a  nitrogen  heteroatom  in  the  conjugated  backbone  of  these 
polymers  presents  a  potentially  severe  breaking  of  both  spatial  and  charge- 
conjugation  symmetry  (CCS),  and  the  addition  of  nonbonding  (n)  orbitals  has 
potentially  major  effects  on  the  photophysics  of  these  systems.  Geometries 
are  optimized  at  the  PM3  Hartree-Fock  level  for  neutral,  singly  charged  and 
doubly  charged  oligomers.  We  find  that  the  geometric  distortions  associated 
with  polaron  formation  are  centered  on  the  vinylene  linkages  in  PPyV-based 
systems  and  on  the  interring  bonds  in  the  PPy-based  systems.  We  discuss 
the  electronic  structure  at  the  PM3  level  applying  configuration  interaction 
between  singly  excited  states  (SCI),  and  we  demonstrate  that  the  lowest- 
lying  (n  — ►  T*)  states  of  the  ideal  polymer  chain  are  well  above  the  lowest 
(x  T*)  states,  predicting  strong  fluorescence  in  these  systems.  Deviations 
from  ideal  geometry,  however,  can  lead  to  substantial  mixing  of  the  (x  — ►  x*) 
and  (n  — ►  x’)  mzinifolds,  thereby  altering  this  conclusion.  We  calculate  ab¬ 
sorption  spectra  for  neutral,  singly  charged  (polaron),  doubly  charged  (bipo- 
laron)  and  triplet-state  oligomers  using  tl>fe  INDO/SCI  technique.  For  PPyV, 
comparison  of  oligomers  with  differing  spatial  symmetry  allows  the  isolation 
of  the  effects  of  CCS  breaking.  All  calculated  spectra  are  in  good  agreement 
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with  eip«imental  results  eud  mdicate  that  the  symmetry  breaking  due  to  the 
nitrogen  heteroalom  is  weak.  In  particular,  the  polaron  mduces  a  two-peak 
in-gap  feature  into  the  absorption  spectrum  and  the  bipolaron  a  single-peak 
feature,  as  is  seen  in  the  analogous  all-hydrocarbon  polymers. 


Typeset  using  REVTg^C 


1.  INTRODUCTION 


With  the  advent  of  the  polymer  light-emitting  device  (LED)  [1,2],  an  understanding  of 
the  electronic  structure  of  conjugated  polymers  has  become  increasingly  important.  The 
system  that  has  received  the  most  attention  is  poly(p-phenylene  vinylene)  (PPV)  and  its 
derivatives,  as  they  are  the  most  promising  candidate  to  date  for  LED  applications  [3]. 
Due  to  the  nature  of  the  electron- phonon  interaction  in  quasi-one-dimensional  systems  [4], 
the  fundamental  charge  excitations  of  PPV  are  the  polaron  (radical  cation/anion)  and  the 
bipolaron  ( iication/dianion).  The  presence  of  such  an  excitation  on  a  ^PV  chain  leads  to 
a  destabilization  of  the  levels  near  the  Hartree-Fock  bandgap,  resulting  in  two  in-gap  levels 
whose  occupation  depends  on  the  charge  of  the  species.  Because  of  the  C2h  spatial  symmetry 
of  PPV,  the  polaron  and  bipolaron  produce  a  two-peak  subgap  feature  and  a  one-peak  subgap 
feature  in  the  absorption  spectrum  for  the  respectively  [5].  It  is  now  generally  accepted  that 
the  fundamental  neutral  excitation  of  the  polymer  is  an  exciton  with  substantial  Coulomb 
binding  energy  [6-10].  The  excitonic  nature  of  neutral  excited  states  requires  the  application 
of  correlated  models  such  as  those  used  in  quantum  chemistry,  and  many  such  studies  have 
been  reported  for  oligomeric  model  systems  for  PPV  [5,9,11-13]. 

While  PPV  possesses  both  charge-conjugation  symmetry  (CCS)  and  Cj/^  spatial  symme¬ 
try,  much  recent  attention  has  been  given  to  systems  in  which  these  symmetries  are  broken. 
For  PPV,  this  can  be  achieved  experimentally  by  the  addition  of  sidegroups  to  the  phenyl 
ring.  Rice  and  Gartstein  have  introduced  a  semi- analytical  model  for  describing  the  ab¬ 
sorption  spectrum  of  these  systems  based  upon  excitation  of  the  benzene  monomer  and  the 
concept  of  correlation  gaps  [14—16].  Their  model  predicts  a  characteristic  four-peak  absorp¬ 
tion  spectrum  for  broken-CCS  systems,  with  the  positions  and  oscillator  strengths  of  the 
four  peaks  determined  by  microscopic  Coulomb  parameters.  A  four-peak  spectrum  is  also 
predicted  by  Chandross  et  al.  [9]  based  on  calculations  within  the  Pariser-Parr-Pople  (PPP) 
model  augmented  with  single-excitation  configuration  interaction  (SCI).  While  the  number 
of  features  predicted  is  independent  of  the  model  considered,  the  predicted  polarization 
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dependence  is  different  for  the  two  models. 

In  this  pnper,  we  consider  the  electronic  structure  of  two  representative  pyridine-based 
polymers,  namely  poly(p-pyridyl  vinylene)  (PPyV)  and  poly{;vpyridine)  (PPy).  shown  in 
Fig.  1.  These  pyridine-based  analogs  of  PPV  and  polyCp-phenylene)  (PPP),  respectively, 
recently  have  shown  promise  for  LED  appHcations  [ir-19|,  and  their  photophysics  shows 
interesting  deviations  from  that  of  PPV  and  PPP  [20,21|.  In  particular,  the  inclusion  of  a 
nitrogen  heteroatom  in  the  unit  cell  of  the  pyridine-based  systems  leads  to  the  introduction 
of  nonbonding  (n)  orbitals  (from  the  N  lone  pair)  into  the  electronic  structure,  in  addition  to 
the  the  breaking  of  CCS  and  spatial  symmetry.  We  assess  the  role  of  these  modificalions  m 
the  electronic  structure  using  various  neglect-of-differential-overlap  (NDO)  methods.  Such 
methods  have  been  applied  with  great  success  to  the  study  of  PPV  oligomers  by  Bredas 
and  co-workers  (5,11-13|.  We  demonstrate  that  many  experimental  observations  on  the 
pyridine-based  systems  can  be  accounted  for  by  our  calculations,  including  the  shape  of 
ground-state  and  photoinduced  absorption  features,  the  strong  fluorescence  properties,  and 
the  morphology  dependence  of  intersystem  crossing.  We  further  demonstrate  that  symmetry 
breaking  has  little  effect  on  the  observed  spectra,  so  that,  e.g.,  the  bipolaron  possesses  a 
single  in-gap  absorption  feature.  We  compare  our  results  with  other  proposed  models  for 

broken-symmetry  systems. 


II.  METHODS 

AU  calculations  were  performed  on  a  Pentium  computer  running  at  60  MHs  using  Hyper- 
Chem,  a  software  package  which  is  commercially  available  from  AutoDesk,  Inc.  The  me 
of  nomenclature  used  for  the  oligomers  considered  in  the  calculations  is  demonstrated  in  Fig. 
1.  OUgomers  are  referred  to  as  ‘PYnXX’  and  ‘PYVnXX,’  where  ‘n’  refers  to  the  number  of 
pyridine  rings  in  the  oUgomer  and  ‘XX’  refers  to  the  regioisomer  considered.  For  the  “head- 
to-tail”  (HT)  regioisomer,  neighboring  pyridine  rings  are  oriented  m  the  same  direction.  For 
the  “head-to-head”  (HH)  regioisomer,  neighboring  rings  are  oriented  oppositely. 
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The  geometries  of  the  oligomers  were  optimized  at  the  PM3  level  [22].  Of  the  three 
neglect-of-diatomic-differential-overlap  methods  (MNDO,  AMI,  PM3),  PM3  gave  the  most 
accurate  geometry  for  the  pyridine  ring  with  respect  to  STO-3G  and  3-21G  ab  initio  calcula¬ 
tions.  In  particular,  AMI  gave  an  overly  distorted  pyridine  ring.  The  neutral  (ground-state) 
and  bipolaron  geometries  were  optimized  at  the  restricted  Hartree-Fock  (RHF)  level  and  the 
polaron  and  triplet  geometries  at  the  unrestricted  Hartree-Fock  (UHF)  level.  The  evolution 
of  the  ground-state  electronic  structure  was  considered  at  the  PM3-I-SCI  RHF  level.  Absorp¬ 
tion  spectra  were  calculated  using  the  INDO/SCI  method  in  the  RHF  or  restricted  open-shell 
Hartree  Fock  (ROHF)  approximation.  The  ‘ZINDO/S’  parameterization  of  INDO  was  used 
[23],  as  is  appropriate  for  spectroscopic  calculations.  The  size  of  the  Cl  space  was  varied 
with  oligomer  size,  as  is  discussed  below.  Experimental  spectra  were  simulated  based  on  the 
calculated  oscillator  strengths,  /,  and  transition  energies,  Eo,  using  a  Gaussian  broadening: 

/  •  exp[— (FI  —  EoY/2cr^],  with  a  =  0.15  eV  (unless  otherwise  indicated). 

III.  RESULTS 

A.  Geometry 

Representative  bond  lengths  for  a  central  repeat  unit  of  PM3-optimized  PYV6HH  and 
PY6HH  are  summarized  in  Fig.  2.  While  no  experimental  measurements  of  these  quantities 
are  currently  available,  the  PM3  results  are  quite  similar  to  those  calculated  for  smaller 
systems  by  ab  initio  techniques,  and  are  close  to  those  reported  for  PPV  oligomers  [12], 
with  a  general  shortening  of  the  bonds  around  the  N  atom.  When  allowed  to  vary,  the 
interring  torsion  angles  between  rings  in  the  PPyV  oligomers  were  found  to  be  very  small  in 
all  cases,  reflecting  the  shallow  steric  potential  in  these  and  related  systems  [24].  Therefore, 
to  save  computational  time,  geometries  for  PPyV  obgomers  were  forced  to  planarity  in 
the  calculations.  For  head-to-head  oligomers  of  PPy,  the  torsion  angle  between  two  facing 
pyridine  rings  (e.g.,  rings  1  and  2  of  PY6HH  above)  was  similarly  found  to  be  essentially  zero 


due  to  lack  of  steric  effects.  Adjacent  “back-to-back”  rings  (e.g.,  rings  2  and  3  of  PY6HH), 
however,  were  found  to  have  a  total  torsion  angle  of  ~  40  degrees  between  them  (i.e.,  ±20 
degrees  as  mesured  from  the  undisorted  plane).  For  head-to-tail  oligomers  of  PPy,  all  rings 
are  found  to  be  separated  by  ~  40  degree  torsion  angles. 

The  geometry  relaxation  phenomena  involved  in  the  formation  of  polaronic  and  triplet 
states  are  demonstrated  for  the  PYV6HH  and  PY6HH  electron  polaron  (P“),  electron  bipo- 
laron  (BP"*)  and  triplet  exciton  (T)  in  Fig.  3,  where  changes  in  bond-lengths  are  shown 
along  a  path  through  the  backbone  of  the  molecule  which  contains  all  of  the  nitrogen  atoms. 
(The  sequence  of  bonds  is  shown  in  inset  of  the  Figure.)  Similar  distortions  are  seen  in 
the  portions  of  the  rings  not  contained  in  this  path.  The  results  for  hole  polarons  and 
bipolarons  are  very  similar  and  are  not  shown.  The  option  of  geometric  relaxation  in  the 
singlet  excited  state  was  not  available  with  the  software  used.  The  geometric  distortions 
are  nearly  independent  of  the  charge  of  the  species,  and  are  strongly  localized  to  the  central 
two  rings  and  their  surrounding  bonds.  For  PYV6HH,  the  vinylene  unit  experiences  bond- 
order  reversal,  with  the  double  bonds  becoming  longer  and  the  single  bonds  shorter,  and 
the  pyridine  ring  becomes  more  quinoid-like.  Notably,  the  polaron  produces  the  weakest 
geometric  distortion,  as  also  is  seen  in  PPV  oligomers  [5].  For  all  species,  the  greatest  geo¬ 
metric  distortion  is  found  to  occur  in  the  central  vinylene  unit,  again  similar  to  conclusions 
drawn  concerning  the  PPV  polaron  [25].  The  vinylene-centered  polaron  explains  the  lack 
of  any  particularly  strong  vibrations  associated  with  the  nitrogen  heteroatom  in  ns  and  ms 
photoinduced  infrared  absorption  experiments  [26] . 

Likewise,  for  PY6HH,  all  distortions  are  centered  around  the  central  interring  bond,  with 
similar  large  distortions  occurring  in  the  two  adjacent  interring  bonds.  The  shortening  of 
these  bonds,  together  with  an  observed  planarization  of  the  molecule,  reflect  aromatic-to- 
quinoid  conversion  upon  polaron  or  triplet  formation  [27].  For  both  oligomers,  the  extent 
of  the  topological  defect  induced  by  the  polaron  or  triplet  state  is  four  to  six  rings,  again 
similar  to  the  results  for  PPV  oligomers  '5l. 

To  demonstrate  that  these  conclusions  are  not  merely  the  result  of  the  symmetry  of 
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the  molecule,  similar  results  are  displayed  for  PYV5HT  and  PY5HT  in  Fig.  4.  Again,  the 
changes  in  bond-length  are  shown  through  a  path  which  includes  all  of  the  nitrogens.  As 
with  the  “hexamers,”  a  vinylene-centered  polaron  is  found  for  the  PPyV  oligomer,  and  the 
di'-tortions  in  the  PPy  oligomer  are  centered  on  the  two  central  interring  bonds.  Both  the 
magnitude  and  extent  of  the  geometric  distortions  are  similar  to  the  hexamer  case,  as  well 
as  to  the  results  found  for  PPV  [5]. 

B.  Electronic  structure  of  neutral  oligomers 

To  tacditate  the  discussion  of  the  electronic  structure  calculations,  the  Huckel  x- 
bandstructures  of  PPyV  and  PPy  are  shown  in  Fig.  5.  This  simplistic  calculation  was 
performed  assuming  a  single  hopping  integral  t  and  a  nitrogen  site  energy  Sn  =  ec  —  0.5t 
[22].  A  head-to-tail  geometry  was  assumed  in  both  cases  (so  that  the  unit  cell  contains 
a  single  pyridine  ring),  and  the  PPy  polymer  was  assumed  to  be  planar.  Nonplanarity 
of  PPy  merely  results  in  a  narrowing  of  the  calculated  Huckel  bands.  The  PPyV  band- 
structure  contains  four  occupied  x  bands  (1-4)  and  four  unoccupied  x*  bands  (1**4*),  as 
shown.  The  PPy  bandstructure  contains  three  occupied  and  three  unoccupied  x— bands. 
The  <r-bandstructures  for  the  polymers  are  not  shown,  but  Ukely  overlap  with  the  outermost 
x-bands  based  on  the  HF  calculations.  A  band  of  states  Tom  the  occupied  nonbonding  (n) 
levels,  is  also  not  shown,  but  presumably  lies  close  to  the  uppermost  occupied  x  band  (1). 
Each  occupied  band  is  approximately  mirrored  by  an  unoccupied  (starred)  band,  indicating 
that  PPyV  and  PPy  are  nearly  charge-conjugation  symmetric  within  this  simplistic  model. 

In  PPV  and  PPP,  the  x-levels  can  be  divided  into  “localized”  and  “delocalized”  orbitals 
[9].  The  localized  orbitals  derive  from  benzene  molecular  orbitals  which  lack  density  at  the 
para  sites  along  which  the  polymer  is  connected.  Therefore,  both  PPP  and  PPV  possess 
one  dispersionless  occupied  band  (1)  and  one  dispersionless  unoccupied  band  (I*).  The 
corresponding  bands  in  PPyV  are  those  labelled  2  and  2*  in  the  Figure;  interaction  with 
the  delocalized  bands  gives  these  bands  a  nonvanishing  dispersion  in  PPyV.  In  PPP,  each 
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localized  band  is  accidentally  degenerate  with  a  delocalized  band  at  fc  =  x/2a;  in  PPy,  these 
accidental  degeneracies  are  broken,  leading  to  well-dispersed  occupied  (1,2)  and  unoccupied 

(1*,  2*)  bands. 

In  the  SCI  calculations  presented  below,  the  Cl  space  was  taken  to  include  transitions 
from  levels  corresponding  to  the  occupied  bands  1,  2,  3  and  n  to  levels  corresponding  to  the 
unoccupied  bands  1*,  2*  and  3*  of  PPyV.  For  PPy  oligomers,  transitions  from  bands  1,  2 
and  n  to  bands  1*  and  2*  were  considered.  At  the  HF  level,  the  tt  bands  4  and  4  of  PPyV 
and  3  and  3*  of  PPy  are  found  to  be  strongly  intermingled  with  cr  bands  and  are  difficult  to 
account  for.  Nevertheless,  calculations  on  PYV2HH  and  PY2HH  which  include  all  of  the  tt 
electrons  suggest  that  the  contributions  of  the  outermost  bands  to  the  Cl  are  minimal. 

We  now  consider  the  evolution  of  the  electronic  structure  with  oligomer  length.  We  first 
consider  the  results  of  SCI  calculations  performed  on  top  of  the  PM3  calculations.  Since  we 
are  considering  only  single  excitations,  the  many-electron  states  calculated  from  the  Cl  can 
be  labeled  as  (o  — ►  where  o  and  u  describe  the  occupied  state  and  unoccupied  state, 

respectively,  from  which  the  transition  is  derived,  and  2S  -f- 1  represents  the  spin  multiplicity 
of  the  final  state  (singlet  or  triplet).  Four  categories  are  allowed  depending  on  the  occupied 
state  from  which  the  electron  is  transferred  and  on  the  spin  multiplicity:  (ir  states 

derive  from  promotion  of  an  electron  from  an  occupied  tt  level  to  an  unoccupied  tt*  level; 

states  derive  from  the  promotion  of  an  electron  from  a  nonbonding  state. 
Note  that  all  nonbonding  states  are  occupied  in  both  polymers.  In  planar  PPyV,  the  four 
categories  of  states  are  rigorously  forbidden  to  mix.  In  PPy,  however,  the  nonplanarity  of  the 
ground-state  molecule  can  lead  to  substantial  mixing  of  (tt  tt*)  and  (n  ->  tt*)  mamfolds. 
Nevertheless,  even  with  40  degree  torsion  angles  between  rings,  the  orbitals  of  PPy  oligomers 
still  largely  retain  their  (n  -  r’)  and  (tt  -  tt*)  characters,  thus  allowing  an  approximate 
categorization  of  the  levels. 

Of  extreme  importance  for  luminescence  in  a  molecule  containing  n  electrons  is  the 
relative  position  of  the  lowest  (n  — *•  ir*)  states  with  respect  to  the  lowest  singlet  (ir  — ♦ 
TT*)  state.  This  is  because  intersystem  crossing  (i.e.,  singlet-to-triplet  conversion)  between 
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(n  —*  TT*)  and  (tt  —>■  x*)  states  is  allowed  to  first  order  in  the  spin-orbit  interaction  [28]. 
Therefore,  molecules  in  which  an  (n  — +  x*)  state  lies  below  the  lowest  (x  — >  x*)^  state 
generally  fluoresce  only  weakly,  if  at  all. 

Figure  6  demonstrates  the  evolution  of  the  lowest  (x  — >  x*)  and  (n  — +  x*)  states  of 
PPyV  and  PPy  oligomers  as  a  function  of  oligomer  length.  Here,  m  refers  to  the  number 
of  non-hydrogen  atoms  in  a  path  connecting  the  two  ends  of  the  oligomer.  Note  that  these 
calculations  were  performed  on  head-to-head  oligomers  in  hopes  of  maximizing  the  interac¬ 
tion  between  the  n  levels.  The  [n  x*)  states  show  very  little  evolution  with  chain  length, 
a  result  which  can  easily  be  rationalized  d  ie  to  the  localized  nature  of  the  n  level.  Also, 
the  singlet-triplet  splitting  is  small  for  these  levels,  as  is  characteristic  of  molecular  systems 
[29].  On  the  other  hand,  the  (x  —*■  x*)  states  exhibit  strong  evolution  with  chain  length. 
The  position  of  the  various  states  in  the  polymer  can  be  estimated  via  the  ansatz  [30] 

^oligomerij^')  —  ^pojj/mer  "i"  ?  (f) 

m 

where  A  is  a  constant.  The  polymer  levels  therefore  correspond  to  the  y-intercepts  in  Fig. 
6.  The  experimentally  determined  location  of  the  lowest  (x  — >  x*)^  state  [26]  is  indicated 
by  the  arrow  in  the  Figure,  and  is  seen  to  be  in  good  agreement  with  the  predicted  value. 
For  PPyV,  it  is  clear  that  the  lowest  (n  x*)  levels  lie  well  above  the  lowest  (x  — >  x*)^ 
level,  therefore  predicting  strong  fluorescence  for  this  system.  For  PPy,  the  (n  — >  x*)  levels 
are  significantly  closer  to,  yet  still  above,  the  lowest  (x  — x*)'  level. 

In  comparing  the  resiilts  of  different-length  oligomers,  it  is  imperative  that  the  Cl  space 
is  kept  size  consistent  [22].  Thus  the  number  of  levels  included  in  the  Cl  calculation  must 
increase  in  proportion  with  the  length  of  the  oligomer.  The  calculated  results  correspond  to 
Cl  spaces  of  size  2N  -i-  1,  where  N—7  x  5  for  the  dimer,  15  x  11  and  23  x  17  for  the  PPyV 
hexamer.  The  full  calculation  for  the  hexamer  was  actually  too  large  for  the  computer  to 
handle,  so  the  Cl  space  was  taken  just  large  enough  to  give  a  proper  calculation  of  the 
(n  — +  X*)  levels  for  this  oligomer.  Note  that  in  a  recent  calculations  for  PPV  oligomers  il2j 
the  position  of  the  lowest  triplet  state  was  found  not  to  evolve  with  chain  length,  suggesting 
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that  the  triplet  is  localized  to  a  single  ring.  In  the  present  (size-consistent)  calculation  for 
PPyV  and  PPy,  this  is  not  the  case:  the  evolution  of  the  lowest  triplet  with  chain  length  is 

almost  as  strong  as  that  of  the  lowest  singlet. 

It  should  be  noted  that  the  results  presented  in  Fig.  6  are  for  the  ideal  (geometry- 
optimized)  polymer.  When  ring-torsional  disorder  is  present,  however,  the  conclusions  may 
be  invahd.  Indeed,  substantial  mixing  of  (ir  ->  tt*)  and  (n  tt*)  levels  is  seen  even  in 
geometry-optimized  PY4HH,  particularly  in  higher-lying  states,  as  is  demonstrated  below. 

A  more  dramatic  example  is  that  of  PY4HH  with  one  of  the  end  rings  twisted  out  of  the 
plane  by  90  degrees.  In  this  case,  the  n-Uke  orbital  of  the  rotated  ring  mixes  in  roughly  equal 
proportion  with  the  x-like  orbital  in  the  more  planar  portion  of  the  molecule.  The  tnplet 
state  formed  between  this  composite  orbital  and  the  LUMO  lies  below  the  lowest  (tt,  tt  ) 
state,  resulting  in  a  potentially  large  intersystem  crossing  rate.  Similar  conclusions  can 
reached  for  PPyV  oligomers.  Experimentally,  a  morphology- dependent  intersystem  crossing 
rate  is  seen,  with  triplet  generation  in  disordered  powder  samples  roughly  ten  times  higher 
than  in  more-ordered  film  samples  or  solutions  [20].  We  propose  the  morphology  dependence 
to  be  due  to  stronger  ring-torsional  disorder  in  powder  samples,  leading  to  enhanced  triplet 
generation  via  (n  — +  x*)  and  (tt  — ►  tt*)  state  mixing. 

C.  Absorption  spectra  of  neutral  oligomers 

We  now  discuss  the  ground-state  absorption  spectrum  of  neutral  PPyV  and  PPy.  For 
these  calculations,  the  INDO/SCI  technique  was  implemented.  It  is  possible  to  calculate 
absorption  spectra  using  the  PM3  technique  Mowed  by  an  SCI  calculation.  While  the 
PM3  HF  electronic  structure  is  probably  a  much  better  description  of  the  actual  electronic 
structure  than  that  obtained  from  INDO,  the  PM3+SCI  spectra  are  generaUy  untrustwor¬ 
thy,  depending  sensitively  on  the  size  of  the  Cl  space.  In  particular,  PM3+SCI  calculations 
for  polaronic  and  triplet  states  proved  to  be  unreliable,  producing  many  negative  eigenval¬ 
ues.  The  INDO/SCI  method,  however,  is  parameterized  specifically  for  the  reproduction  of 
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spectroscopic  results  and  is  much  more  reliable.  Unfortunately,  due  to  the  relatively  poor 
description  of  the  HF  electronic  structure,  the  oscillator  strengths  of  low-lying  transitions 
are  somewhat  exaggerated,  and  the  entire  spectrum  is  in  general  blue-shifted  from  the  ex¬ 
perimental  result.  Previous  authors  have  corrected  for  these  shortcomings  by  correcting  the 
intensities  based  on  valence-effective  Hamiltonian  calculations  and  by  scaling  the  spectrum 
along  the  energy  axis  [5].  We  choose  to  merely  present  the  “raw”  INDO/SCI  result  here. 

While  the  main  goal  of  the  following  is  to  reproduce  experimental  observations,  we  also 
attempt  to  asses  the  role  of  both  spatial  and  CCS  breaking  in  optical  spectra.  In  general, 
violatio"'  of  either  symmetry  will  lead  to  shifting  of  oscillator  strength  '  and  the  appearance 
of  new  transitions  in  optical  spectra.  Ideally,  one  could  assess  the  role  of  CCS  breaking 
by  examining  the  Cl  coefficients  associated  with  various  transitions.  For  example,  if  PPyV 
were  charge-conjugation  symmetric,  then  transitions  between  Hiickel  bands  1  and  2*  of  Fig.  5 
would  have  corresponding  degenerate  transitions  between  bands  2  and  1*.  Therefore,  any  Cl 
expansion  which  contains  contributions  from  a  1  — >  2*  transition  should  necessarily  contain 
an  equal  contribution  from  a  2  ^  1*  transition,  and  the  extent  to  which  one  transition  is 
favored  over  the  other  in  the  Cl  expansion  should  be  a  manifestation  of  CCS  breaking. 

Unfortunately,  it  is  a  general  fzuling  of  NDO  techniques  that  the  HF  electronic  structure 
is  not  charge-conjugation  symmetric  even  in  charge-conjugation  symmetric  systems  such  as 
PPV  oligomers.  This  conclusion  is  evidenced  by  our  own  results  on  PPV  and  PPP  oligomers, 
which  reproduce  the  results  of  Cornil  et  al.  [11].  Similar  phenomena  have  been  ascribed  to 
the  poor  description  of  the  unoccupied  levels  at  the  INDO  level  [5].  While  it  is  possible 
that  this  failing  may  be  compensated  by  the  Cl  calculation,  it  is  nevertheless  impossible  to 
directly  assess  the  role  of  CCS  breaking  in  the  pyridine-based  systems  by  simply  examining 
Cl  coefficients. 

On  the  other  hand,  the  nature  of  the  PPyV  system  makes  it  possible  to  separate  the  two 
symmetries  by  considering  the  two  regioisomers  independently.  This  is  because  the  head-to- 
head  oligomers  retain  the  Can  spatial  symmetry  of  PPV  while  demonstrating  broken  CCS. 
On  the  other  hand,  the  head-to-tail  ohgomers  do  not  possess  either  symmetry.  One  can 
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therefore  evaluate  the  role  of  CCS  breaking  by  comparing  head-to-head  spectra  with  those 
of  PPV  oligomers,  and  the  effects  of  additional  spatial  symmetry  breaking  can  be  assessed 

by  comparing  head-to-head  and  head-to-tail  spectra. 

Figure  7  shows  the  calculated  absorption  spectrum  for  PYV4HH  (upper)  and  PYV4HT 
(lower).  The  two  spectra  are  quite  similar,  indicating  that  the  spatial  symmetry  breabng 
introduced  by  the  nitrogen  is  minimal.  The  soUd  and  dashed  spectra  correspond  to  Gaussian 
broadening  factors  of  0.15  eV  and  0.35  eV,  respectively.  The  0.35  eV-broadened  spectra  best 
represent  the  types  of  Unewidths  seen  in  experiment.  These  spectra  consist  of  four  peaks, 
labeled  I-IV  in  the  Figure,  together  with  a  high  energy  feature  (HE).  The  four-peak  spectrum 
is  seen  in  experiment  [31],  with  the  HE  peak  out  of  experimental  range,  and  the  intensities 
and  positions  of  the  peaks  are  reproduced  weU  by  the  calculation  to  within  a  scaling  factor 
in  energy.  In  the  spectrum,  as  in  experiment,  peaks  I  and  IV  are  most  prominent,  with  peak 

III  merely  a  shoulder  to  peak  IV. 

In  comparison  with  calculated  results  for  PPV  tetramers  [11],  the  PYV4HH  spectrum 
contains  several  extra  transitions  due  to  broken  CCS.  A  four-peak  absorption  spectrum  is 
predicted  for  polymers  with  broken  CCS  in  the  model  of  Gartstein  et  al.  [15]  and  within 
the  PPP  model  calculations  of  Chandross  et  d.  [9].  In  the  present  calculation,  the  presence 
of  four  peaks  in  the  spectrum  is  intimately  .elated  to  the  choice  of  broadening  parameter 
used.  Indeed,  the  spectrum  under  0.15  eV  broadening  appears  to  show  at  least  6  peaks, 
in  addition  to  the  HE  feature.  One  difficulty  arises  from  finite-size  effects:  Some  of  the 
transitions  represented  in  the  tetramer  spectrum  will  disappear  in  the  polymer  spectrum 
due  to  momentum  conservation  in  the  translationaUy  invariant  system.  If  spectra  of  various 
length  oUgomers  are  compared,  it  is  found  that  the  transitions  in  the  bands  labeled  I,  III 
and  IV  are  robust-their  relative  oscillator  strengths  and  positions  are  relatively  independent 
of  oligomer  length.  The  number  of  transitions  in  the  band  labeled  II,  however,  is  strongly 
dependent  on  oUgomer  length,  indicating  the  role  of  finite  size  effects. 

One  obvious  difference  between  the  two  regioisomer  spectra  is  that  the  head-to-tail  spec¬ 
trum  contains  many  more  significant  transitions  than  the  head-to-head  spectrum.  As  a 
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result  of  the  sum  rule  [32],  the  total  oscillator  strength  in  the  two  spectra  should  be  the 
same.  Therefore,  the  average  oscillator  strength  of  the  transitions  in  the  head-to-tail  isomer 
is  somewhat  less  than  that  in  the  head-to-head  isomer.  This  is  a  manifestation  of  spatial 
symmetry  breaking  due  to  the  nitrogen:  the  head-to-head  isome’-  is  inversion  symmetric, 
whereas  the  head-to-tail  isomer  is  not.  Therefore,  there  are  fewer  forbidden  transitions  in  the 
head-to-tad  spectrum.  Nevertheless,  the  general  shape  of  the  two  spectra  are  very  similar, 
suggesting  spatial  symmetry  breaking  plays  only  a  weak  role  in  the  spectrum. 

We  now  consider  in  detail  the  origin  of  the  features  in  the  PYV4HH  spectrum.  It  is 
useful  to  consider  the  transitions  in  terms  of  the  Htickel  bandstructure  presented  above. 
Comparison  of  the  wavefunctions  of  the  HF  molecular  orbitals  with  those  expected  from 
the  Huckel  bandstructure  allows  direct  association  of  the  molecular  orbitals  of  the  tetramer 
with  the  bands  1-4  and  i*-4*  of  Fig.  5.  Likewise,  the  nonbonding  orbitals  of  the  N  atoms 
can  be  separated  out  from  the  tt  and  tt*  levels.  For  PYV4HH,  the  first  three  unoccupied 
HF  levels  correspond  to  band  1*  of  the  polymer,  the  next  four  to  2*,  and  the  following  four 
to  3*.  Note  that  there  may  be  some  ambiguity  in  this  assignment  due  to  the  fact  that  the 
“tetramer”  is  two  carbons  short  of  four  complete  unit  cells;  however,  the  HF  orbitals  show 
a  clear  distinction  between  the  first  three  unoccupied  orbitals  and  the  next  four  unoccupied 
orbitals.  A  similar  level  scheme  exists  for  the  occupied  orbitals,  with  the  four  nonbonding 
orbitals  interjected  periodically.  Within  the  Huckel  theory,  one  expects  the  lowest  lying 
transition  to  result  from  the  1  — >  1*  transition.  Two  additional  low-energy  features  are 
expected  due  to  the  1  — ►  2*  and  2  — 1*  transitions. 

Table  I  shows  the  squared  Cl  coefficients  for  the  various  transitions  in  the  calculated 
absorption  spectrum.  For  ease  of  presentation,  the  individuaJ  levels  involved  in  the  transition 
are  not  singled  out;  only  the  bands  between  which  the  transitions  occur  are  indicated.  For 
example,  the  peak  labeled  ‘I’  in  the  PYV4HH  spectrum  is  largely  (81%)  due  to  the  HOMO 
—>■  LUMO  transition.  This  result  is  typical  of  suc!h  calculations  i  ll].  Of  the  remaining  19%, 
16%  of  the  oscillator  strength  of  this  transition  comes  from  other  1  — v  1*  transitions;  e.g., 
HOMO-1  -+  LUMO-i-1.  The  feature  labeled  TP  contains  a  large  contribution  from  1  — »  2* 
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Md  2  1-,  as  expected.  Note  that  substantial  mixing  with  the  1  ^  V  transition  occurs  in 

tUs  peak.  Peaks  ‘IIP  and  ‘IV’  demonstrate  strong  mixing  between  various  bands.  The  HE 
feature  originates  largely,  though  not  entirely,  from  Iransitidhs  Involving  the  3  and  3’  bands. 
If  these  bands  are  left  out  of  the  calculation,  the  HE  band  disappears  from  the  spectrum, 
leaving  the  other  bands  unchanged.  Interestingly,  configurations  involving  bands  3  and  3’ 


also  play  a  substantial  role  in  peaks  III  and  IV. 

Both  the  model  of  Gartstein  el  al.  and  the  calculations  of  Chandross  el  ol.  predict  strong 

mixing  between  the  various  Hiickeblike  bands;  therefore,  a  direct  comparison  between  these 
models  and  the  present  calculations  based  on  Cl  coefficients  is  quite  difficult.  The  two 
models  differ  substantiaUy,  however,  in  their  predictions  of  the  polarization  dependence  of 
the  various  features  in  the  absorption  spectrum.  In  particular,  ihe  model  of  Gartstein  el  ai. 


predicts  that  peak  II  is  polarized  perpendicular  to  (but  in  the  plane  of)  the  chain;  whereas 
Chandross  el  al.  predict  polarization  along  the  chain.  For  this  reason,  we  have  tabulated  the 
polarization  dependence  of  the  observed  features  in  Table  I.  Much  of  the  osciHator  strength 
of  the  peak  II  is  polarized  along  the  chain,  as  predicted  by  the  PPP  calculation.  Caution 
should  be  taken  in  this  interpretation,  however,  as  finite-size  effects  may  also  play  a  role  in 


the  polarization  dependence  [16]. 

Finally,  before  leaving  the  PYV4HH  spectrum,  we  note  that  the  nonbonding  orbitals 
participate  only  weakly  in  the  observed  features.  The  lowest-lying  (n  -e  W)  transition  occurs 
at  4.583  eV,  and  its  osciUator  strength  is  less  than  one  percent  of  that  of  the  first  optical 
transition.  Note  that  its  polarization  is  labeled  ‘i,’  indicating  polarization  perpendicular  to 
the  plane  of  the  molecule.  No  other  (n  ^  x’)  state  contributes  significantly  to  the  spectrum. 

Figure  8  displays  the  INDO/SCI  calculated  absorption  spectra  tor  the  two  regioisomers 
of  PY4.  Again,  the  sobd  line  corresponds  to  a  broadening  of  0.15  eV  and  the  dashed  li 
0.35  eV.  The  0.35  eV  spectrum  best  reproduces  the  experimental  Unewidths.  This  spectrum 
shows  two  clear  peaks,  labeled  1  and  II  in  the  Figure,  and  a  high-energy  feature,  HE.  Again, 
the  predicted  two-peak  spectrum  closely  reproduces  the  experimental  result  :3l!,  with  the 
HE  peak  out  of  the  range  of  the  experiment. 
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In  PY4HH  there  are  two  pairs  of  coplanar  rings  with  a  torsion  angle  of  40  degrees 
between  them.  In  PY4HT,  all  rings  have  a  40  degree  torsion  angle  between  them.  The 
resulting  change  in  delocalization  leads  to  a  blueshift  in  the  low-energy  feature  in  going 
from  the  head-to-head  to  the  head-to-tail  spectrum.  Again  the  two  spectra  are  quite  similar 
suggesting  the  weakness  of  spatial  symmetry  breaking.  As  in  the  case  of  PYV4,  the  PY4HT 
spectrum  displays  a  greater  number  of  transitions  than  the  PY4HH  spectrum.  Note  that, 
however,  PY4HH  is  nof  inversion  symmetric  due  to  the  40  degree  torsion  angle  in  the  center  of 
the  molecule.  It  is  possible  that  the  PY4HH  molecule  is  approximately  inversion  symmetric, 
however,  ?gain  leading  to  a  larger  number  of  forbidden  transitions  th?"'  in  PY4HT. 

Table  II  shows  the  squared  Cl  coefficients  for  PY4HH.  Again,  we  label  the  levels  according 
to  their  association  with  the  Hiickel  bands  of  Fig.  5.  As  with  PYV4,  peak  I  in  the  PY4HH 
spectrum  is  almost  entirely  1  — ►  1*.  The  other  features  are  again  strong  admixtures  of  the 
various  transitions.  The  features  in  the  intermediate  region  between  peaks  I  and  II  contain 
strong  contributions  from  1  — ^  1*,  1  — ►  2*,  and  2  — ►  1*,  again  as  expected  from  Hiickel 
theory.  The  feature  labeled  II  contains  contributions  from  2^2*  and  the  high-energy  side 
of  the  1  — >  1*  transition.  The  features  I,  II  and  HE  are  robust  as  the  length  of  the  oligomer 
is  changed.  The  number  and  size  of  the  features  contributing  to  the  region  between  I  and 
II  is  dependent  on  oligomer  length,  indicating  finite-size  ettects. 

It  is  interesting  to  note  the  role  of  the  (n  — ►  tt*)  state  in  the  nonplanar  PY4HH  molecule. 
One  effect  of  the  nonplanarity  is  to  mix  the  (n  tt*)  and  (tt  — ►  tt*)  states.  Therefore,  many 
of  the  transitions  in  Table  II  contain  contributions  from  both  types  of  states.  As  a  result  of 
this  mixing,  the  (n  x*)  states  acquire  substantial  oscillator  strength.  This  is  dramatically 
demonstrated  by  the  first  two  transitions  identified  with  peak  II.  Both  of  these  transitions 
sire  roughly  50%  (n  — ►  x*)  and  possess  oscillator  strengths  near  unity.  Furthermore,  the 
polarization  of  the  (n  — ►  x*)  state  is  now  in  the  ^plane’  of  the  molecule  (defined  as  passing 
through  the  long  axis  of  the  molecule  and  lying  halfway  between  the  two  planes  define  by  the 
pairs  of  rings  in  PY4HH).  On  the  other  hand,  some  (x  x*)  transitions  are  now  polarized 
perpendicular  to  the  plane  of  the  molecule  (again  denoted  by  ‘z’  in  the  Table). 
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D.  Polarons,  bipolarons  and  triplet  excitons 


We  now  turn  to  a  discussion  of  the  absorption  spectra  of  the  polaron,  bipolaron  and 
triplet  exciton.  At  the  HF  level,  the  introduction  of  a  polaron  or  exciton  to  a  neutral  chain 
results  in  a  destabilization  of  the  HOMO  and  LUMO  levels,  bringing  them  into  the  HF 
bandgap.  This  is  represented  schematicaUy  in  Fig.  9.  The  two  levels  I  and  u  are  split  of 
from  the  bands  1  and  1*,  respectively.  The  Figure  demonstrates  the  case  of  an  electron 
polaron.  For  the  case  of  the  electron  bipolaron,  the  u  level  is  doubly  occupied.  For  the 
case  of  the  triplet  exciton,  both  levels  are  singly  occupied.  It  is  interesting  to  note  that 
destabiliL.  tion  of  the  levels  around  the  bandgap  occurs  at  the  HF  .cvd  even  before  the 
molecule  is  allowed  to  relax  geometricaUy.  The  geometric  distortion  from  polaron  formation 
further  destabilizes  these  levels. 

The  spectra  for  the  electron  polaron  and  bipolaron  on  PYV4HH  and  PYV4HT  are  shown 
in  Fig.  10.  The  spectra  of  the  two  regioisomers  are  very  similar,  again  indicating  that  spatial 
symmetry  breaking  is  weak  in  the  pyridine-based  systems.  The  squared  Cl  coefficients  are 
Usted  for  the  head-to-head  oHgomer  in  Tables  III  and  IV  for  the  polaron  and  bipolaron, 
respectively.  From  Hiickel  theory,  one  expects  up  to  three  optical  transitions  for  the  electron 
polaron  in  charge-conjugation  symmetric  systems,  namely  u  ^  1*,  /  ^  u,  and  1  ^  1*.  A 
fourth  transition,  1  u,  is  degenerate  with  the  /  1*  transition  in  charge-conjugation 

symmetric  systems,  but  can  result  in  an  extra  feature  in  non-charge-conjugation  symmetric 
systems.  The  bipolaron  can  have  up  to  two  in-gap  transitions  (u  ^  1*,  /  ^  1*),  independent 
of  CCS.  In  PPV  oligomers,  the  /  ^  1*  and  1  u  transitions  are  forbidden  by  the  Cih  spatial 
symmetry,  leading  to  a  two-peak  polaron  and  one- peak  bipolaron  spectrum  [5]. 

The  two-peak  polaron  and  one-peak  bipolaron  spectral  signatures  are  also  seen  in  the 
PPyV  tetramers.  Unlike  in  the  case  of  PPV,  however,  in  the  PPyV  oligomers  the  high-energy 
polaron  feature  and  the  single  bipolaron  feature  are  actually  composed  of  several  transitions. 
These  transitions  are  indicated  in  Tables  III  and  IV.  In  the  head-to-head  spectrum,  PI  is 
largely  (61%)  due  to  transitions  from  the  upper  polaron  level  to  the  new  LUMO  (u  ^  1*), 
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as  expected  [33].  Of  the  two  largest  transitions  which  contribute  to  P2  (1,586  eV  and 
1.728  eV),  only  one  is  the  expected  I  u  transition.  The  second  transition  is  actually  the 
u  2*  transition  (1.728  eV).  As  PYV4HH  possesses  the  same  C2h  spatial  symmetry  as 
PPV  oligomers,  we  suggest  that  this  extra  transition  results  from  the  breaking  of  charge- 
conjugation  symmetry.  The  /  — ►  1^  and  1  — >  u  transitions  which  are  forbidden  in  PPV 
oligomers  [5]  also  do  not  occur  in  the  PYV4HH  spectra,  consistent  with  the  preservation  of 
C2h  spatial  symmetry  in  this  oligomer.  The  absence  of  a  third  feature  in  the  head-to-tail 
spectrum  indicates  that  spatial  symmetry  breaking  is  weak.  While  it  is  also  true  that  the 
/  — ►  1*  transition  is  predicted  to  be  weak  in  systems  where  polaron  confinement  is  strong 
[33],  we  note  that  the  two-peak  polaron  spectrum  can  be  seen  in  PPyV  oligomers  even  before 
geometric  relaxation  is  performed.  Therefore,  we  believe  that  the  absence  of  the  third  peak 
is  truly  due  to  symmetry  considerations,  and  not  due  to  polaron  confinement. 

In  the  bipolaron  spectrum,  the  singular  feature  in  the  spectra  is  due  almost  entirely 
(90%)  to  transitions  from  the  upper  bipolaron  level  to  the  LUMO  {u  — ►  1*),  as  expected 
[33] .  Again,  the  I  1*  transitions  which  are  forbidden  in  PPV  oligomers  [5]  do  not  occur 
in  PPyV  oligomers.  In  particular,  a  second  feature  is  absent  from  the  head-to-tail  spec¬ 
trum,  again  suggesting  the  weak  role  of  spatial  symmetry  breaking.  Again,  the  single-peak 
bipolaron  spectrum  can  be  recovered  in  PPyV  oligomers  before  geometric  relaxation  has  oc¬ 
curred,  indicating  that  the  absence  of  a  second  peak  is  related  to  the  relatively  weak  spatial 
symmetry  breaking  in  these  systems,  and  not  due  to  polaron  confinement  considerations. 

We  also  calculated  spectra  for  hole  polarons  and  bipolaxons.  The  shape  of  these  spectra 
are  similar  to  those  of  their  electron  counterparts  (i.e.,  two-peak  polaron  spectrum,  etc.), 
with  the  spectrum  of  the  positively  charged  entity  rigidly  redshifted  from  its  negatively 
charged  counterpart  by  several  tenths  of  an  eV .  While  it  is  tempting  to  ascribe  this  red- 
shift  to  the  broken  charge  conjugation  symmetry  in  the  system,  our  calculations  on  PPV 
oligomers  indicate  that  electron  and  hole  polarons  have  different  spectra  even  in  this  charge- 
conjugation  symmetric  system.  Again,  this  result  can  be  explained  by  the  poor  description 
of  the  unoccupied  states  afforded  by  the  INDO  calculation  i5]. 
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The  polaron  and  bipolaron  spectra  are  shown  for  the  two  regioisomers  of  PPy  in  Fig. 
11.  The  squared  Cl  coefficients  are  presented  in  Tables  V  and  VI,  respectively.  The  two- 
peak  polaron  spectrum  and  one-peak  bipolaron  spectrum  are  reminiscent  of  PPP  oligomers 
[34],  indicating  only  weak  symmetry  breaking.  The  polaron  spectra  are  somewhat  more 
greatly  affected  by  regioregularity;  however,  the  bipolaron  spectra  are  quite  similar.  For 
the  head-to-head  spectra,  again  the  PI  transition  is  largely  (80%)  upper-polaron-level  to 
LUMO,  as  expected.  In  contrast  to  the  case  of  PYV4HH,  in  the  PY4HH  spectrum  P2 
contains  essentiaUy  one  transition  (2.817  eV)  that  is  72%  /  ^  u,  as  expected.  Again,  the 
main  bipolaron  transition  (BPl)  is  almost  entirely  (92%)  u  ^  1*,  as  expected.  Unlike  in  the 
case  of  PYV4,  in  PY4  the  bipolaron  possesses  a  substantial  shoulder  (BP2).  This  shoulder 
is  mainly  due  to  transitions  to  the  2*  level.  Again,  the  /  1*  and  1  «  transitions  that 

are  forbidden  in  PPV  oUgomers  do  not  contribute  to  the  spectra  of  the  PPy  oligomers. 

Experimentally,  a  two-peak  feature  is  seen  in  the  millisecond  photoinduced  absorption 
spectrum  of  both  PPyV  and  PPy,  particularly  in  dims  [31].  From  the  discussion  above, 
these  spectra  are  assigned  to  polaron  absorptions;  however,  exactly  which  types  of  polarons 
(electron  or  hole,  head-to-head  or  head-to-tail)  are  present  is  not  determinable  from  the 
above  considerations.  Bipolarons  have  not  been  clearly  detected  in  any  of  the  pyridine- 
based  polymers. 

The  triplet-triplet  spectra  for  PYV4  and  PY4  are  shown  in  Figs.  12  and  13,  respectively. 
The  squared  Cl  coefficients  for  the  head-to-head  oUgomers  are  shown  in  Tables  VII  and 
VIII  respectively.  Again,  the  spectra  are  fairly  independent  of  the  regioisomer  considered, 
suggesting  a  weak  role  of  spatial  symmetry  breaking.  For  PPV,  it  has  been  indicated  that 
a  single  strong  transition  dominates  the  triplet -triplet  spectrum  [12].  Our  own  calculations 
also  indicate  a  single  strong  transition  in  oUgomers  of  PPP.  For  the  triplet  m  PYV4  and 
PY4,  as  with  the  polaron  and  bipolaron,  the  main  feature  is  composed  of  several  transitions, 
each  of  which  is  a  strong  admixture  of  various  HP  configurations.  The  additional  transitions 
in  the  PYV4HH  spectrum  again  reflect  the  breaking  of  charge- conjugation  symmetry  m  this 
system.  .As  the  triplet  spectrum  is  Ukely  more  affected  by  correlations  than  the  polaron  or 
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bipolaron  spectra,  it  is  difficult  to  compare  the  observed  transitions  with  those  expected 
from  Huckel  theory. 

Experimentally,  a  single  feature  is  seen  [31]  whose  position  and  shape  is  matched  quite 
w^Il  by  the  PYV4HT  spectrum  when  the  results  are  scaled  with  the  absorption  spectra.  On 
the  other  hand,  the  PPy  triplet-triplet  spectrum  is  observed  experimentally  to  be  nearly 
identical  to  the  PPyV  spectrum.  This  observation  is  not  reproduced  by  the  calculation. 

IV.  SUMMARY  AND  CONCLUSIONS 

The  role  of  spatial  symmetry  breaking  in  the  pyridine- based  systems  studied  is  at  best 
weak.  This  is  evidenced  by  a  general  lack  of  dependence  of  calculated  spectroscopic  results 
on  the  regioisomer  considered.  In  particular,  comparison  of  the  two  regioisomers  of  PPyV 
allows  the  direct  evaluation  of  spatial  symmetry,  as  the  head-to-head  oligomers  retain  the  C2h 
spatial  symmetry  of  PPV  oligomers;  whereas  the  head-to-tail  oligomers  show  broken  spatial 
symmetry.  For  the  PPyV  oligomers,  spatial  symmetry  breaking  leads  to  additional  high- 
energy,  low-oscillator  strength  transitions  in  the  absorption  spectrum  of  neutral  oligomers; 
however,  the  general  shape  of  the  spectrum  is  independent  of  regioisomer.  Polaron,  bipolaron 
and  triplet  spectra  are  somewhat  less  affected  by  spatial  symmetry  breaking,  which  results 
in  slight  shifts  in  oscillator  strength  zunong  the  various  transitions. 

The  role  of  CCS  breaking  is  somewhat  more  elusive,  as  the  INDO  technique  generally 
gives  inadequate  results  even  for  charge- conjugation  symmetric  systems.  Nevertheless,  the 
effects  of  CCS  breaking  in  the  pyridine- based  systems  can  be  isolated  by  comparing  the 
results  for  head-to-head  oligomers  of  PPyV  with  those  of  PPV,  as  both  possess  the  same 
C2/1  spatial  symmetry.  Although  the  breaking  of  CCS  introduces  several  new  transitions  in 
the  pyridine- based  systems,  particularly  in  the  absorption  spectrum,  in  general  the  spectro¬ 
scopic  properties  of  the  these  systems  closely  resemble  those  of  their  hydrocarbon  analogues, 
indicating  that  CCS  breaking  is  also  weak.  In  particular,  the  electron  polaron  introduces 
two  in-gap  features  into  the  absorption  spectrum  of  the  pyridine- ba.sed  oligomer,  and  the 
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electron  bipolaron  introduces  a  single  in-gap  feature,  as  predicted  for  PPV  [5].  The  general 
shapes  of  the  calculated  absorption  spectra  for  the  pyridine- based  oUgomers  are  consistent 

with  that  expected  from  theories  in  which  CCS  is  weakly  broken  [16]. 

The  (n,7r*)  states  play  a  unique  role  in  these  systems.  The  present  calculations  find  these 
levels  too  high  in  energy  to  affect  the  photophysics  of  the  ideal  geometry-optimized  polymer, 
predicting  strong  luminescence  for  these  systems.  However,  nonplanarity  can  induce  mixing 
of  these  levels  with  the  (7r,7r*)  states,  allowing  them  to  lower  their  energy  and  increase 
their  osciUator  strength  substantially.  We  therefore  predict  a  morphology  dependence  to 
the  photophysics,  in  particular  to  rates  of  intersystem  crossing. 

Indeed,  the  morphology  dependence  of  intersystem  crossing  is  seen  m  experiment  [31|. 
Several  other  experimental  observations  are  predicted  well  by  the  current  calculations,  no¬ 
tably,  the  four-peak  PPyV  absorption  spectrum,  the  two-peak  PPy  absorption  spectrum, 
and  the  two-peak  polaron  absorption  spectra.  In  addition,  the  PPy  V  triplet-triplet  spectrum 
is  well  predicted  by  our  calculations.  Finally,  the  vinylene-centered  polaron  predicted  by  the 
calculations  is  consistent  with  photoinduced  infrared  absorption  data,  which  demonstrate  a 
lack  of  ring-based  modes. 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research. 
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FIGURES 

FIG.  1.  Repeat  units  for  PPyV  and  PPy;  representative  oUgomers  used  in  this  study.  Nomen- 
clature  is  described  in  the  text. 

FIG.  2.  Bond  lengths  (in  Angstroms)  of  central  repeat  unit  of  neutral  PYV6HH  (upper)  and 
PY6HH  (lower)  after  PM3  geometry  optimization. 

FIG.  3.  Ch».g.  in  bond  length  vt.  bond  number  for  PYV6HH  (upper)  and  PY6HH  (lower). 
The  geometrical  distortion!  are  mapped  along  a  path  containing  all  N  atoms,  as  is  indicated  in  the 

insets. 

FIG.  4.  Change  in  bond  length  vs.  bond  number  for  PYV5HT  (upper)  and  PY5HT  (lower). 
The  geometrical  distortions  are  mapped  along  a  path  containing  aU  N  atoms,  as  is  indicated  in  the 

insets. 

FIG.  5.  Calculated  Hiickel  x-electron  structure  for  PPyV  (upper)  and  PPy  (lower). 

FIG.  6.  Evolution  of  energies  of  lowest  (ir  -  :r‘)  and  (n  t*)  states,  calculated  at  the 
PM3+SCI  level,  as  a  function  of  inverse  oligomer  length  (in  number  of  carbons)  for  PPyV  (upper) 
and  PPy  (lower)  head-to-head  oUgomers.  Arrows  represent  experimentally  determined  {x  x*)i 

location. 

FIG.  7.  INDO/SCI  calculated  absorption  spectrum  of  neutral  PYV4HH  (upper)  and  PYV4HT 
(lower),  assuming  a  0.15  eV  Gaussian  broadening  (soUd  Unes)  or  0.35  eV  Gaussian  broademng 

(dashed  lines). 

FIG.  8.  INDO/SCI  calculated  absorption  spectrum  of  neutral  PY4HH  (upper)  and  PY4HT 
(lower),  assuming  a  0.15  eV  Gaussian  broademng  (soUd  Unes)  or  0.35  eV  Gaussian  broademng 

(dashed  lines). 

FIG.  9.  Schematic  representation  of  Hiickel  levels  of  an  electron  polaron  on  a  PPyV  chain. 
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FIG.  10.  INDO/SCI  calculated  absorption  spectrum  of  an  electron  polaron  (solid)  and  electron 
bipolaron  (dashed)  on  PYV4HH  (upper)  and  PYV4HT  (lower). 

FIG.  11.  INDO/SCI  calculated  absorption  spectrum  of  an  electron  polairon  (solid)  and  electron 
bipolaron  (dashed)  on  PY4HH  (upper)  tind  PY4HT  (lower). 

FIG.  12.  INDO/SCI  czJculated  triplet- triplet  absorption  spectrum  of  PYV4HH  (upper)  tind 
PYV4HT  (lower). 

FIG.  13.  INDO/SCI  calculated  triplet-triplet  absorption  spectrum  of  PY4HH  (upper)  and 


PY4HT  (lower). 


TABLES 

TABLE  I.  INDO/SCI  transitions  energies,  polarizations,  oscillator  strengths,  and  squared  Cl 
coefficients  for  neutral  PYV4HH.  Polarization  represents  the  angle  of  the  transition  dipole  (in 
deg.)  in  the  plane  of  the  molecule  with  respect  to  the  long  axis  of  the  molecule. 

TABLE  n.  INDO/SCI  transitions  energies,  polarizations,  oscillator  strengths,  and  squared  Cl 
coefBcients  for  neutral  PY4HH. 

TABLE  in.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients 
for  an  electron  polaron  on  PYV4HH. 

TABLE  rv.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients 
for  an  electron  bipolaron  on  PYV4HH. 

TABLE  V.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients  for 
an  electron  polaron  on  PY4HH. 

TABLE  VI.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients 
for  an  electron  bipolaron  on  PY4HH. 

TABLE  Vn.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients 
for  triplet  eiciton  on  PYV4HH. 

TABLE  Vni.  INDO/SCI  transitions  energies,  oscillator  strengths,  and  squared  Cl  coefficients 
for  triplet  eiciton  on  PY4HH. 
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Figure  1,  Blatchford,  Gustafson  and  Epstein 
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Figure  5,  Blatchford,  Gustafson  and  Epsteii 
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Figure  6,  Blatchford,  Gustafson  and  Epstein 
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Figure  7,  Blatchford,  Gustafson  and  Epstein 
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Figure  8,  Blatchford,  Gustafson  and  Epstei 


Figure  9,  Blatchford,  Gustafson  and  Epstein 
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Figure  10,  Blatchford,  Gustafson  and  Epstein 
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Figure  11,  Blatchford,  Gustafson  and  Epstein 
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Fiqure  12,  Blatchford,  Gustafson  and  Epsteii 
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Figure  13,  Blatchford,  Gustafson  and  Epstein 


i-eak 


III 


IV 


Transition 
energy 
(eV)  _ 


3.496 


4.513 

4.583 

4.640 

5.186 

5.365 


6.016 

6.111 

6.163 

6.271 

6.310 


Pol. 


6.789 

6.864 

6.890 

7.193 

7.248 

7.359 

7.390 

7.497 

7.552 

7.706 

7.862 

7.958 

8.135 

8.210 


16.9 


18.7 

z 

71.1 

17.3 

13.9 


10.3 

12.3 

57.1 

60.2 

71.2 


63.1 

43.3 

55.9 

47.3 

53.3 

15.0 

11.9 
39.5 

48.9 


Oscillator 

Strength 


2.7514 


0.2709 

0.0152 

0.0808 

0.3319 

0.2150 


0.0889 

0.0642 

0.0667 

0.0154 

0.4419 


85.0 

39.2 
25.4 

32.2 


0.0886 

0.5761 

0.3695 

0.3989 

0.2100 

0.8747 

0.0419 

0.0815 

0.0506 


Squared  Cl  Expansion  Coefficients 


0.97  (1-^1*) 


0.76  (1^2*) +  0.12  (2^1*) 

0.74  (n^l*)  +  0.25  (n^2*) 

0.65  (l^r)  +0.19  (1-+2*)  +  0.06  (2->1*)  +  0.05  (2->2*) 
0.52  (1  -+2*)  +  0.29  (1  ->r)  +  0.09  (2->1  *) 

0.79  (1->r)  +  0.09  (1-»2*) 


0.1466 

0.7591 

0.0771 

0.3322 


0.43  (1->2*)  +  0.27  (1^1*)  +  0.22  (2-+1*) 

0.47  (1-+3*)  +  0.28  (1-+r)  +  0.16  (1-^2*) 

0.33  (1^1*)  +  0.28  (1^3*)  +  0.23  (2^1*)  +  0.09  (1-+2*) 

0.52  (2-^1*)  +0.16  (2^1*)  +  0.14  (1-+r)  +  0.07  (1->3*)  + 
0.05  (2-^2*) 

0.40  (2^1*)  +  0.35  (1-+2*)  +  0.0^  (1^1*)  +  0-06  (2-+2p - 

0.23  (3->r)  +  0.20  (1^2*)  +  0.20  (1-^3*)  +0.13  (2-^2*)  + 

0.12  (2->1*)  +0.05  (1^1*) 

0.35  (2-+1  *)  +  0.25  (1^3*)  +  0.20  (1  -^2*)  +  0.07  (2-+2*)  + 
0.06  (1  -+1  *) 

0.36  (1-»3*)  +0.16  (2-»2*)  +0.12  (3^1*)  +  0.11  (1-^2*)  + 
0.11  (2-+r) +  0.09  (1-+1*) 

0.45  (1^2*)  +  0.22  (2^1*)  +  0.12  (2->2*)  +  0.10  (1-+3*) 

0.34  (2^1*)  +  0.23  (2-^2*)  +  0.16  (1-+3*)+  0.10  (1-+2’)  + 
0.06  (3^r) +0.06  (1-+1*) 

0.31  (2-^2*)  +  0.21  (2->r)  +  0.19  (3->1*)  +  0.09  (1->3  )  + 
0.07  (1  ^2*)  +  0.05  (3->2*) 

0.29  (1  -+3*)  +  0.22  (1  ^2*)  +  0.1 3  (1  ^1  *)  +  0.1 2  (3-+1  *)  + 
0.10  (2^2*)  +  0.07  (2-+1*) 

0.28  (1->3*)  +  0.21  (1-^2*)  +  0.15  (2-+2*)  +  0.09  (2-+r)  + 
0.09  (1  -+r)  +  0.06  (2^3*)  +  0.05  (3^2*) 

0.44  (2-^2*)  +  0.11  (2-+3*)  +  0.10  (1-»1*)  +  0.10  (1-+3  )  + 

0.09  (3-+2*)  +  0.06  (3^1*)  +  0.06  (2^1*) _ _ 

0.28  (1-+1*)  +  0.18  (1-+3*)  +  0.15  (3->r)  +  0.12  (1-^2  )  + 
0.10  (2->r)+ 0.07  (2-+2*) 

0.37  (2-»r)  +  0.16  (3-^r)  +  0.15  (1^2*)  +  0.11  (2-+3  )  + 
0.11  (2->2*) 

0.26  (2-+r)  +  0.25  (1-^2*)  +0.15  (2-+2*)  +  0.14  (2-^3  )  + 
0.07  {3->2’)  +  0.06  (1  ^3*) 

0.31  (1->3*) +  0.23(3^1*)  + 0.19  (2^3*) +  0.11  (2-^r)  + 
0.05  {2-+2*) 

0.23  (1-+3’)  +  0.18  (3-^r)  +  0.16  {2-+2*)  +  0.08  (3-+3  )  + 

0.10  (2->3*)  +  0.08  (2-+r)  +  0.08  (3->2*) 

+  0.06(1-+2*) 
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8.309 

8.450 

8.572 

8.602 

8.620 

8.741 

8.766 

8.872 

8.910 

9.093 

9.242 

9.332 

9.688 

9.696 


63.9 

60.1 

79.0 

83.2 

38.7 

64.2 

55.1 

45.1 

71.8 

62.3 

61.1 

62.9 

58.1 

65.2 


0.5434 

0.0173 

0.0162 

0.8939 

0.0769 

0.1182 

0.0898 

0.0835 

0.1017 

0.0720 

0.0775 

0.0325 

0.0420 

0.0324 


0.31  (2^2*)  +  0.24  (2^1*)  +  0.23  (3^2*)  +  0.06  (3-»r) 
0.34  (1-^2*)  +  0.28  (2-^2*)  +  0.1 1  (2-^1*)  +  0.09  (3^2*)  + 
0.07  (1-^r) 

0.28  (2^2*)  +  0.26  (2^1*)  +  0.24  (3^1*)  +  0.14  (2^3*)  + 
0.08  (3-»2*) 

0.21  (1^3*)  +0.18  (3^1*)  +  0.17  (2^2*)  +  0.16  (3^3*)  + 
0.14  (2^3*)  +  0.08  (2->1*)  +  0.06  (3^2*) 

0.43  (2^2*)  +  0.24  (1^3*)  +  0.12  (2^3*)  +  0.07  (1-»2*)  + 
0.05  (3^1*)  +  0.05(3-^2*) 

0.36  (2-^2*)  +  0.25  (2->3*)  +  0.16  (2^1*)  +  0.06  (1^3*)  + 
0.05  (3^2*) 

0.39  (2->2*)  +  0.35  (3-^2*)  +  0.10  (1-^3*)  +  0.06  (3->1*) 
0.32  (3-+2*)  +  0.26  (2^3*)  +  0.17  (2^2*)  +  0.12  (2-^1*)  + 
0.06  (3-^1*) 

0.36  (1  ^3*)  +  0.18  (3^2*)  +  0.14  {2->r)  +  0.14  (3^1*)  + 
0.12  (3-^3*) 

0.35  (1^3*)  +  0.34  (2^3*)  +  0.10  (2-+r)  +  0.06  (3^3*)  + 
0.06  (2-+2*)  +  0.05  (3-+2*) 

0.37  (3^2*)  +  0.17  (2^3*)  +  0.10  (1^3*)  +  0.1 1  (2^2*)  + 
0.10  (2^1*) +  0.05  (3->3*) 

0.29  (2^2*)  +  0.28  (3-^1*)  +  0.20  (2^3*)  +  0.06  (2->r)  + 
0.05  (3->3*)  +  0.05  (1->3') 

0.61  (2->2*)  +  0.16  (3-^2*)  +  0.08  (2-^3*)  +  0.07  (3^1*)  + 
0.05  (2-+r) 

0.26  (2^3*)  +  0.18  (3^3*)  +  0.18(2^2*)  +  0.16  (3^1*)  + 
0.09  (3^2*)  +  0.05  (2^1*) 


Table  I  cont’d,  Blatchford,  Gustafson  and  Epstein 


Peak 

Transition 

energy 

(eV) 

Pol. 

Oscillator 

Strength 

Squared  Cl  Expansion  Coefficients 

1 

i _ \  - 

4.100 

3.8 

1.7212 

0.95(1^1*) 

4.447 

z 

0.0129 

0.84  (n^r)  +  0.07  (n-^2*)  +  0.05  (1^1*) 

4.486 

14.9 

0.0154 

0.83  (n^r)  +  0.08  (n  >2*)  +  0.06  (1-^1*) 

4.783 

24.9 

0.1294 

0.77  (1->r)  +  0.10  (2^1*)  +  0.06  (2^2*) 

5.026 

42.7 

0.1838 

0.49  ( 1  ^2*)  +  0.32  ( 1  1  *)  +  0. 1 1  (2^  1  *)  +  0.06  (2-^2*) 

5.599 

29.4 

0.0130 

0.48  (n^r)  +  0.44  (n^2*) 

5.660 

z 

0.0200 

0.76  (1-41*)  +  0.09  (n-^r)  +  0.08  (1-^2*) 

5.930 

35.4 

0.2712 

0.60  (1-^1*)  +  0.29  (1^2*)  +  0.05  (2^1*)  +  0.05  (n->1*) 

6.310 

42.9 

0.0337 

0.43  (1^2*)  +  0.36  (W1*)  +  0.13  (2^1*) 

6.690 

61.1 

0.1188 

0.42  (1  ^1  *)  +  0.29  (1-42*)  +  0.21  (2^1  *) 

6.701 

z 

0.0543 

0.34  (1->2*)  +  0.33  (Wr)  +  0.26  (2^1*) 

6.730 

84.8 

0.2833 

0.51  (1^2’)  +  0.28  (2-41*)  +  0.13  (1^1*) 

6.804 

65.1 

0.1573 

0.51  (1  ^1*) +  0.30  (2^1*) +  0.07  (1^2*)+ 0.06  (n-^1*) 

II 

7.091 

16.2 

1.1261 

0.39  (n^l*)  +  0.19  (2-41*)  +  0.16  (1^1*)  +  0.15  (2->2*)  + 

0.10  (n^2*) 

7.115 

13.1 

0.8210 

0.29  (2^1*)  +  0.27  (n^r)  +  0.15  (n-^2*)  +  0.09  (2-^2*)  + 

0.06  (1^2*) +  0.05  (1-^1*) 

7.144 

z 

0.0119 

0.47  (n^r)  +  0.34  (n-^2*)  +  0.08  (1-42*)  +  0.04  (2-^1*) 

7.210 

z 

0.0897 

0.37  (1-^2*) +  0.33  (2^1*) +0.11  (1^1*) +  0.11  (n^l*) 

7.231 

64.4 

0.8140 

0.36  (2^1*)  +  0.32  (1^2*)  +  0.12  (2-42*)  +  0.12  (I-4I*) 

7.443 

68.3 

0.2265 

0.29  (n-4l*)  +  0.24  (n^2*)  +  0.16  (2-41*)  +  0.13  (1-^1*)  + 

0.07  (2-42*) +  0.09  (1-^2*) 

7.515 

69.3 

0.3541 

0.34  (2^1*)  +  0.34  (I-4I*)  +  0.15  (W2*)  +  0.08  (n->2*) 

7.589 

4.4 

0.0376 

0.29  (1^1*)  +  0.20  (n->r)  +0.18  (n->2*)  +  0.17  (2-41*)  + 
0.10(2^2*) 

7.783 

79.4 

0.0484 

0.27  (n-42*)  +  0.27  (W2*)  +  0.13  (n-»1*)  +  0.12  (1^1*)  + 

0.10  (2-^2*)  + 0.09  (2->1*) 

HE 

7.948 

46.6 

0.0959 

0.28  (2-^1*)  +  0.21  (n-4l*)  +  0.20  (1^2*)  +  0.13  (n-42*)  + 

0.09  (2-42*)  +  0.06  (1^1*) 

7.988 

0.9 

0.1152 

0.45  (2-41*)  +  0.42  (2-42*)  +  0.06  (1  ->1*) 

8.040 

68.6 

0.3116 

0.41  (n-42*)  +  0.25  (2^1*)  +  0.14  (n-4l*)  +  0.12  (1-^2*) 

8.226 

64.2 

0.1460 

0.28  (1-42*)  +  0.22  (2-42*)  +  0.16  (2^1*)  +  0.14  (1^1*)  + 
0.10  (n-4l*)  + 0.07  (n-42*) 

8.331 

67.6 

0.0936 

0.29  (1^2*)  +  0.27  (2-41*)  +  0.19  (I-4I*)  +  0.18  (2^2  ) 

8.342 

84.2 

0.7729 

0.35  (2-41*)  +  0.34  (2-42*)  +0.20  (1-42*)  +0.11  (I-4I*) 

8.353 

z 

0.0743 

0.37  (2^1  *)  +  0.29  (2-42*)  +  0.21  (1  -42*)  +  0.1 0  (1 -4I  *) 

8.427 

66.6 

0.0180 

0.64  (n-4r)  +  0.28  (n-42*)  ^  «  -to  m  i*\ 

8.479 

z 

0.0238 

0.27  (1-42*)  +  0.24  (2-42*)  +  0.20  (n^l  )  +  0.12  (I-4I  )  + 
0.11  (2^1*)  +  0.06  (n-42*) 

8.514 

74.8 

0.4391 

0.39  (2^1*)  +  0.37  (2^2*)  +  0.19  (1^2*) 

8.747 

87.5 

0.3627 

0.40  (2-42*)  +  0.24  (2-41*)  +  0.21  (I-4I*)  +  0.12  (W2*) 

8.784 

60.1 

0.1131 

0.41  (1^2*)  +  0.26  (2-42*)  +  0.25  (2-41*)  +  0.06  (I-4I*) 

8.962 

34.7 

0.0455 

0.30  (2^2*)  +0.29  (2^1*)  +0.20  (n^2*)  +  0.11  (n-4l*)  + 
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8.981 

43.1 

0.2060 

9.532 

74.2 

0.0355 

9.711 

86.3 

0.0304 

0.07(1^2*) 

0.46  (n^2*)  +  0.20  (n^l*)  +  0.14  (2^2*)  +  0.12  (2-^1*)  + 
0.05(1^2*) 

0.48  (2^2*)  +  0.29  (2-^1*)  +  0.17  (1-^2*) 

0.55  (1-^2*)  +  0.34  (2^2*)  +  0.06  {n-^2*) 


Table  II  cont’d,  Blatchford,  Gustafson  and  Epstein 


T  ransition 
energy 
(eV) 

Oscillator 

Strength 

Squared  Cl  Expansion  Coefficients 

0.871 

0.7266 

0.61  (u-^r)  +  0.16  (u^2*)  +  0.1 1  (3-^2*)  +  0.05  (u-^S*) 

1.281 

0.0625 

0.89  (u — >2  ) 

1.586 

0.8837 

0.24  (^-^u)  +  0.22  (u^r)  +  0.17  (u-»2*)  +  0.15  (3-+2*) 

+  0.05  (^-^1*) 

1.728 

0.6858 

0.54  (u-+2*)  +  0.19  (^-+u)  +  0.07  (3-»2*) 

2.283 

0.0380 

0.43  (u->3*)  +  0.15  (^-+•2*)  +  0.10  (3^3*)  +  0.07  (1-+3*) 

2.783 

0.1975 

0.37  (^-^u)  +  0.36  (3->2*)  +  0.06  (u-+3*) 

2.871 

0.0557 

0.39  (3-43*)  +  0.15  (^-42*)  +  0.10  (3^1*)  +  0.08  (u^3*) 

+  0.06  (^->u) 

Table  III,  Blatchford,  Gustafson  and  Epstein 


T  ransition 
energy 
(eV) 

Oscillator 

Strength 

Squ'irGd  Cl  Expansion  Coofficionts 

_ iZ-Li - 

1.552 

4.0342 

0.90  (u->1*)  +0.06  (^-4l*) 

2.624 

0.0318 

0.80  (u -42*)  +0.12  (^-42*) 

2.794 

0.1549 

0.62  (u-43*)  +  0.28  (u^2*)  +  0.06  (^-42*) 

Table  IV,  Blatchford,  Gustafson  and  Epstein 


Tranjition 

energy 

..(eV) 

Oscillator 

Strength 

Squared  Cl  Expansion  Coefficients 

0.837 

0.4631 

0.88  (u->1*) 

1.639 

0.0807 

0.92  (u->2*) 

2.817 

1.1926 

0.72  (^^u)  + 0.06  (1^1*) 

Table  V,  Hlatchford,  Gustafson  and  Epstein 


T  ransition 
energy 
(eV) 

Oscillator 

Strength 

Squared  Cl  Expansion  Coefficients 

1.814 

2.8016 

0.92  (u-»1*) 

2.253 

0.0179 

0.84  (u-^r)  +  0.08  (u-»2*) 

2.451 

0.4837 

0.86  (u-»2*) +  0.12  (u^r) 

Table  VI,  Blatchford,  Gustafson  and  Epstein 


T  ransition 
energy 
(eV) 

Oscillator 

Strength 

Squared  Cl  Expansion  Coefficients 

0.726 

0.0579 

0.40  (u^r)  +  0.14  (1-^^)  +  0.11  (^-^1*)  +  0.08  (1-4U)+ 
0.07  (u-^3*) 

1.235 

0.0150 

0.79  (u^2*)  +  0.07  ii-^2*) 

1.835 

0.6606 

0.18  +  0.24  (1-^r)  +  0.11  (u-^1*)  +  0.10  (u-^3*) 

+  0.07  (n->1*) 

2.169 

1.0450 

0.25  (u->1*)  +  0.15  (^-+1*)  +  0.14  {i-^3*)  +  0.11  (u-+3*) 

+  0.05  (^-+2*) 

2.434 

0.2195 

0.39  (1— >4  +  J-11  ('“^1*)  +  0.08  (u— >2*)  +  0.05  (1— >2  ) 

2.508 

0.3045 

0.32  (u-^2*)  +  0.19  (1^4  +  0.15  (1-^2*)  +  0.12  (n^2*) 

+  0.07  (^-»2*) 

2.860 

0.1075 

0.19  (u^3*)  +  0.10  (^-+r)  +  0.09  (1-+u)  +  0.09  (1-^1*) 

+  0.06  (u-+2*)  +  0.05  (n-+^)  +  0.05  (3-^^)  +  0.05  (n^l*) 

+  0.05  (1^2’)  +  0.05  (1-^3*) 

3.027 

0.1909 

0.18  (1->u)  +  0.14  (u->3*)  +  0.07  (<;-+7*)  +  0.13  (1->3*) 

+  0.06  (3^1*)  +  0.05  (1  ^1  *)  +  0.05  (n^2*) 

Table  VII,  Blatchford,  Gustafson  and  Epstein 


Transition 

energy 

(eV) 

Oscillator 

Strength 

j  Squared  Cl  Expansion  Coefficients 

1.141 

0.0325 

0.60  (u^r)  +  0.24  (1-»^)  +  0.05  (u-^2*) 

1.959 

0.0702 

0.59  (u-^2*)  +  0.23  (u^l*)  +  0.07  (l^^-) 

2.101 

'  0.0226 

0.50  (u^r)  +  0.37  (u^2*)  +  0.05  (^-^1*) 

2.717 

0.6636 

0.33  (u^r)  +0.32  (U-+2*)  +  0.18  (^-+1*)  +  0.05  (n^^) 

2.906 

1.2306 

0.39  (u^2*)  +  0.35  (1->^)  +  0.14  (u-+r) 

3.330 

0.0601 

0.49  (u^r)  +  0.11  (u^2*)  +  0.09  {^-+1*^  +  0.07  (1^1*) 

+  0.06(1 -.2*) 

3.501 

0.6221 

0.51  (2^^)  +  0.11  (^^1*)  +  0.05  (u-»2*)  +  0.05  (1-»r) 

Table  Vlli,  Biatchford,  Gustafson  and  Epstein 


